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Abstract 

A simple procedure is developed to simultaneously eliminate multiple scattering contributions 
from two reflection electron energy loss spectra (REELS) measured at different energies or for 
different experimental geometrical configurations. The procedure provides the differential inverse 
inelastic mean free path (DIIMFP) and the differential surface excitation probability (DSEP). 
The only required input parameters are the differential cross section for elastic scattering and 
a reasonable estimate for the inelastic mean free path (IMFP). No prior information on surface 
excitations is required for the deconvolution. The retrieved DIIMFP and DSEP can be used to 
determine the dielectric function of a solid by fitting the DSEP and DIIMFP to theory. Eventually, 
the optical data can be used to calculate the (differential and total) inelastic mean free path and 
the surface excitation probability. The procedure is applied to Fe, Co and Ni and the retrieved 
optical data as well as the inelastic mean free paths and surface excitation parameters derived from 
it are compared to values reported earlier in the literature. In all cases, reasonable agreement is 
found between the present data and the earlier results, supporting the validity of the procedure. 
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1. INTRODUCTION. 



The response of a solid to an external electromagnetic perturbation is decribed by the 
dielectric function e(u,q), where fkv and hq are the energy and momentum transfered dur- 
ing the interaction. Knowledge of the dielectric function of a solid is important for many 
branches of physics. The dielectric function can be measured by probing a solid surface 
with elementary particles, e.g. by photons J5| or electrons [6] . With the advent of 

density functional theory beyond the ground state jj], ah initio theoretical calculations of 
optical data have recently become available Q, 9), 10 1. 

From the experimental point of view, reflection electon energy loss experiments represent 
a particularly attractive method to probe the dielectric response of a solid, since the experi- 
ment is very simple. However, quantitative interpretation of the experimental results is not 
straightforward since inside the solid, the electrons experience intensive multiple scattering 
implying that experimental data need to be deconvoluted before information on the dielec- 
tric response of the solid can be extracted from them. In the past, the algorithm of Tougaard 
and Chorkendorff (hereafter designated by " TC" ) [Hi ] has been frequently used for this pur- 



pose. However, it has recently been shown 12J, [13) that the loss distribution retrieved by the 



TC-algorithm is not unique, i.e., it constitutes a not very well defined mixture of so-called 
surface and volume electronic excitations of any scattering order, which makes quantitative 
interpretation and extraction of the dielectric function troublesome. An algorithm which 
is free of these deficiencies has been proposed recently by the present author (denoted as 



"W" -algorithm hereafter) and was succesfully applied to several materials |8|, O, [12|, Il3 |. 
The W-algorithm is more involved than the TC-algorithm and needs evaluation of many 
terms for which the expansion coefficients are tedious and time consuming to obtain. In the 
present work, a simplification of the W-algorithm, (denoted by "SW" -algorithm) is devel- 
oped, following a suggestion made earlier in this connection by Vicanek [14| and succesfully 
tested using REELS spectra of Fe, Co and Ni. The resulting optical data and quantities 
derived from them such as the mean free path for inelastic electron scattering are in good 
agreement with data found in the literature, supporting the validity of the procedure. 
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2. DECONVOLUTION OF REELS SPECTRA. 



A REELS spectrum is made up of electrons that have experienced surface (designated by 
the subscript "s" in the following) and bulk (subscript "b") excitations a certain number of 
times. The number of electrons reaching the detector after participating in n s surface and 
rib bulk excitations is given by the partial intensities A nbUa . Since the occurence of surface 
excitations is localized to a depth region smaller than, or comparable to, the elastic mean 
free path, the par tial intensities for surface and bulk scattering are uncorrelated to a good 



approximation 



15|: 



A nb n a — A nb X A na (1) 

The bulk partial intensities can be obtained most conveniently by means of a Monte Carlo 



calculation, by calculating t 
solid and using the formula 



re distribution of pathlengths Q{s) the electrons travel in the 



16|: 



A nb = / Q(s)W nb (s)ds (2) 



where W n (s) is the stochastic process for multiple scattering: 

w.{.) = (f )»2ELfAi) (3) 

and Xi is the inelastic mean free path. The only physical quantity needed for the calculation 
of the pathlength distribution is the differential cross section for elastic scattering 17| . which 



can be calculated ab initio for free atoms if solid state effects are weak for the considered 
application. The pathlength distribution for a reflection geometry is always much broader 
than the distribution W nb (s) for any value of n^, implying that the reduced partial intensities 
a nb = A nb /A nb=0 depend only very weakly on the value of the inelastic mean free path. 

The surface partial intensities are also assumed to be governed by Poisson statistics. 
Since an electron's trajectory through the surface scattering zone is rectilinear to a good 
approximation (i.e. the pathlength distribution strongly resembles a delta-function), the 
reduced partial intensities for surface scattering a ns = A n JA ns=0 are given by the simple 
equation 

«n s = (n s ) n °/n s \ (4) 
where (n s ) is the average number of surface excitations taking place during reflection, i.e. 



the incoming and outgoing surface crossing is combined in Eqn. 



Commonly, an 



expression of the form: 

is adopted for average number of surface excitations [18j where a s is the so-called surface 
excitation parameter, E is the electron energy and \ii = cos 9i and fi Q = cos8 are the polar 
directions of surface crossing for the incoming and outgoing beam. 

Multiplying the number of electrons arriving in the detector a nbris with the energy dis- 
tribution after experiencing (n&, n s ) collisions, the partial loss distributions, r„ b „ s (T) = 



w 



''^(T — T") (g> w^(T), and summing over all scattering orders, the energy loss spectrum 



oo oo 



y{T) is obtained as [121 ]: 

2/ CO = E E ^^(T - T') ® w^\T) (6) 

n s =0 ri6=0 

Where T denotes the energy loss, the symbol "(g)" represents a convolution over the energy 
loss variable and the quantities and wi are the (n& — l)-fold and (n s — l)-fold self- 
convolution of the normalized energy loss distribution in a single bulk and surface excitation 
respectively, the so-called differential inverse inelastic mean free path (DIIMFP, Wb(T)) and 
differential surface excitation probability (DSEP, w s (T)). Note that the elastic peak needs 
to be removed from a REELS spectrum before analysis, implying that «o,o — [l2j |. 

The subject of the present paper is the retrieval of the quantities u>&CO and w s (T) 
from experimental REELS spectra and determination of the dielectric function from these 
quantities. Recalling the convolution theorem, it is immediately obvious that in Fourier 
space, the spectrum is given by a bivariate power series in the variables Wb(T) and w s (T). 
This implies that a unique solution of these quantities cannot be found by reverting the 
series Eqn. ([6]) since a single equation with two unknowns has no unique solution. However, 
when two loss spectra yi(T) and y% (T) with a different sequence of partial intensities ot nhrLs 
and ni . n , are measured, reversion of the bivariate power series becomes possible using the 
formulae QQ: 

oo oo 

w b{T) = £^«U P) (T-T')®#(T') 

p=0 q=0 

oo oo 

= T,Ty P ,y?(T-T')®y^(T'), (7) 



p=0 q=0 



The above expression constitutes the W-algorit 



and u s „ „ can be obtained as outlined in Refs. 



im for which the required coefficients u 



12 



b 

p,g 



13|- 



The W-a 
materials 12 



gorithm has been succesfully applied to REELS data of a large number of 
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191 ] and the dielectric 



unction of several materials was succesfully extracted 
from the resulting, DIIMFP and DSEP \^,^. However, this approach may be improved upon 
in two respects: first of all, convergence of the series Eqn. (J7|) is relatively slow, implying that 
many terms need to be calculated; secondly, calculation of the higher order coefficients u PA 
is tedious and time consuming. The moderate convergence behaviour of the W-algorithm is 
illustrated in Figure [T^, that shows the various stages of the deconvolution. The first collision 
order n indicated for the different curves implies that the second collision order has been 
taken into account up to n — 1, so, for example, the curve labelled (3,0) is a superposition 
of the (1, 0), (0, 1), (1, 1), (2, 0), (2, 1), (2, 2) and (3, 0)-th order cross-convolutions of the two 
experimental spectra. It is seen that more than 6 scattering orders need to be taken into 
account for the W-algorithm to converge over a loss range of about 100 eV. 

The improvements addressed above can be effected following a suggestion by Vicanek[E 
who pointed out that the TC-algorithm (for the univariate case) can in fact be regarded as 
a lowest order rational fraction expansion (a so-called Pade approximation) of the power 
series expansion of the spectrum. Extension of the concept of the Pade approximation to 



20|. 



multivariate power series, such as Eqn. (J9]) is straightforward 

The simplified algorithm (designated by "SW" -algorithm hereafter) is obtained (in 
Fourier space, indicated by the ("~" )-symbol) by making the rational fraction ansatz : 

N N 

E E a Kiv\y\ 

~ k = 0l = /Q\ 

w = — IT ( } 

E E hjVivi 

fc=CH=0 

where N is the order of the approximation and with &o,o — 1 and ao.o = 0. Note that 
the same formula holds for the quantities w s (T) and Wb{T), only the expansion coefficients 
are different. Explicit expressions for the first order Pade coefficients as well as the system 
of linear equations determining the second order expansion coefficients, are given in the 
appendix. Multiplying Eqn. ([HD with the denominator of the right hand side, and going 
back to real space one finds the deconvolution formula to be given by a Volterra integral 
equation of the second kind (being most harmless numerically): 

N N T r N N 

w ( T ) = J2J2 a ^ Yk ^- / EE 6 *.'^ 1 ^ '-T')w{T')dT' (9) 

fc=0 1=0 T'=0 k=0 1=0 



where the quantity Ykj(T) is the (k, l)-th order cross convolution of the two REELS spectra: 

Y Kl {T) = yf\T-T')®yf{T') (10) 

Since w(T = 0) = the integration on the right hand side of Eqn. (Q can always be carried 
out over the energy loss range for which the loss distribution is already known. For more 



details on the (trivial' 



kind see, e.g., Ref. (211 ] . 



numerical treatment of the Volterra integral equation of the second 



The performance of the SW-algorithm is illustrated in Figure [Tb, that shows the 8-th 
order polynomial expansion approximation (Eqn. <^j) as open circles. The solid curves la- 
belled (1) and (2) are the first and second order Pade approximation results. It is seen 
that the latter converges as good as the polynomial expansion over the first ~150 eV. Even 
the first order approximation leads to reasonable agreement. However, at least the second 
order mixed term (one surface and one bulk excitation) is important to attain quantitative 
agreement. The reason for the good performance of a rational fraction approximation com- 
pared to the polynomial expansion is clear: for a polynomial expansion, for any value of the 
variables a certain power in the series will always be the dominating term and may become 
arbitrarily large. For a rational fraction expansion, a large value of the dominant term in 
the enumerator can always be balanced by an appropriately chosen coefficient of the power 
in question in the denominator. 

The second order SW-algorithm was applied REELS spectra of Ag, Al, Au, Be, Bi, C, 
Co, Fe, Ge, Mo, Mn, Ni, Pb, Pd, Pt, Si, Ta, Te, Ti, V, W and Zn) and was found to 
converge over the measured energy loss range of 140 eV. It can therefore be concluded that 
for an energy loss range of this order the second order SW-algorithm is sufficiently accurate. 
The main merit of the first order approximation is that the expansion coefficients can be 
given in a tractable analytical form (see the Appendix), providing detailed insight in the 
physical behaviour of the deconvolution. For example, it is seen in Eqn. ( TT9]) and (J2TJ) 
that the retrieved DIIMFP is completely independent of the value of the surface excitation 
parameter a s used for the deconvolution. On the other hand, the expansion coefficients 
for the surface loss probability are seen to scale linearly with a s . This can be seen by 
substituting Eqn. (CQ) in Eqn. ( fl9l) and (J2TJ) and using Eqn. (]!]) and Eqn. (]5]). Thus, no 
information on surface excitations whatsoever is needed as input to the procedure (except 
the functional from of the dependence of the surface excitation probability on the energy 
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and surface crossing direction, such as Eqn. (jSJ)). Furthermore, as pointed out above, the 
value of the inelastic mean free path will mainly affect the absolute value of the partial 
intensities, A nb , while the reduced partial intensities a nb change negligibly when the IMFP 
is varied by up to 30%. Reasonably accurate knowledge of the shape and magnitude of 
the elastic scattering cross section is important, however, since the shape of the pathlength 
distribution depends on it. Near a deep minimum of the cross section this dependence may 



even be critical 



16j . giving rise to qualitatively different sequences of partial intensities. 



The higher order expansion coefficients also exhibit the scaling properties addressed 
above. Therefore, using a set of input parameters as discussed above, the normalized DI- 
IMFP is correctly returned by the procedure in absolute units, while the correct shape of 
the DSEP is also obtained, but the absolute value of the DSEP scales with the value of a s 
used in the procedure. 

Finally, it should be pointed out that by a suitable choice of the partial intensities (in 
particular letting (3 nbns — > 0, mimicking the use of a single REELS spectrum) Eqn. (jUJ) 
reduces to the algorithms proposed by Vicanek [yj] or Tougaard and Chorkendorff jll ]. 



3. RETRIEVAL OF OPTICAL DATA FROM THE SINGLE SCATTERING DIS- 
TRIBUTIONS. 



The procedure to extract the dielectric function of a solid from the DSEP and the DIIMFP 
has been outlined in Ref. 9|: The theoretical expressions for the DIIMFP and DSEP are 
fitted to the corresponding experimental results using a suitable model for the dielectric 
function. The differential inverse inelastic mean free path is related to the dielectric function 
e(o>, q) of the solid via the well known formula 221 ] : 

?+ 



W h (u) 



1 

7r£ 



Im 



E[UJ 



L -i dq 
q)f q' 



(11) 



For parabolic bands, the momentum transfer q confined by g_ and q + given by: 



q± 



2E± \/2(E - uj) 



(12) 



Note that atomic units are used in this section. The differential inverse mean free path 
normalized to unity area is related to the unnormalized DIIMFP via Wb{T) = \iWt,(T). 
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For the DSEP, the expression by Tung and coworkers [18|] is used: 

W s (u, 9, E) = P+(u, 9, E) + P~(u, 9, E), (13) 

where the quantity Pt{oJ, 9, E) is defined as 

i+ 

P ° { ^ E) ~ J ie(u J ,q)(e(u,q) + l)l ' (M) 

9- 

and 



( TtttH cosfl±( — ==-)sm0. (15) 

J V Z.C/ 



/2£ 

The normalized differential surface excitation probability is calculated via w s (T,9) = 
W s (T,9)/(n s (9,E)). 

The electromagnetic response of the solid is described in terms of the fit-parameters fi, uji 
and 7j using a model for the dielectric function in terms of a set of Drude-Lindhard oscillators 
for the evaluation of the theoretical expressions for the DIIMFP and DSEP |23j: 



{u 2 - uJi(q) 2 ) 2 + u 2 ^ 2 



ea(w,g) = V ^ f'T^ 2 . — ^ (16) 

^ (u; 2 - uji{q) 2 y + uj 2 ^f 

where u = T is the energy loss and q is the momentum transfer. The static dielectric 
constant is given by e&, /j represents the oscillator strength, ji the damping coefficient and 
Ui is the energy of the 2-th oscillator. Those values of the Drude-Lorentz parameters in 
the above equation that minimize the deviation between the experimental data and the 
theoretical expressions for the DIIMFP and DSEP are taken to parametrize the dielectric 
function. A quadratic dispersion LOi{q) = uii + q 2 /2 was used for the transition described 
by the i— th oscillator, except for those oscillators with an energy higher than the most 
loosely bound core electrons, that are assumed to be described by localized states without 
dispersion. Note that a dielectric function of the form Eqn. (fT6l) implicitly satisfies the 
Kramers-Kronig dispersion relationships. 

The (normalized) distributions Wb(T) and w s (T) obtained from the experimental data 
are simultaneously fitted to the above theoretical expressions by minimizing the function: 

x bXDiiMFp(Pb, fi, ^, 7i) + X s x 2 DS Ep(Ps, fi, W<, 7*)- ( 17 ) 



S 



Here the function \ 2 is the least squares difference between theory and experiment, defined 
in the usual way and X& and X s are weight factors chosen as X& = 1.0 and X s = 0.1, empha- 
sizing the DIIMFP in the fitting procedure. The fit-parameters p^ and p s are multiplicative 
scaling factors for the DIIMFP and the DSEP that compensate an eventual mismatch of the 



surface excitation parameter and IMFP used as input. Using the TPP-2M formula [24j to 
estimate the IMFP, the value of pb deviated from unity by less than 5% for all cases studied, 
while the returned value of p s scales with the value of the surface excitation parameter a s , 
exactly as anticipated in the previous section. 

4. EXPERIMENTAL 

The procedure used to acquire the experimental data has been described in detail before 



251 ]. REELS data were taken for polycrystalline Fe, Co and Ni surfaces in the energy 
range between 300 and 3400 eV for normal incidence and an off-normal emission angle of 
60°, using a hemispherical analyzer operated in the constant-analyzer-energy mode giving 
a width of the elastic peak of 0.7 eV. Count rates in the elastic peaks were kept well below 
the saturation count rate of the channeltrons and a dead time correction was applied to 
the data. The sample was cleaned by means of sputtering with 3 keV Ar + ions. Sample 
cleanliness was monitored with Auger electron spectroscopy. 

For each material the optimum energy combination for the retrieval procedure of two loss 
spectra was determined by inspection of the partial intensities, choosing those energies for 
which on the one hand the determinant A in Eqn. (TSUI) is reasonably large, while, on the 
other hand, those energies were avoided for which the scattering geometry corresponds to 
a deep minimum in the elastic cross section. For these cases it is more difficult to obtain 
the realistic shape of the pathlength distribution since the elastic cross sections are not 
accurately known for such scattering angles and the true electron optical detector solid 
angle also plays a significant role there. In Figure [21 the experimental spectra used in the 
present work are shown as noisy curves. For all selected energy combinations, the shape of 
the loss spectra is seen to be quite similar, but in all cases, a significant difference in the 
relative contribution of surface and bulk excitations is seen, as evidenced by the difference 
in spectral shape below ~20 eV. 

Removal of the elastic peak was achieved by fitting the elastic peak to a combination of a 



9 



Gaussian and a Lorentzian peakshape. Subsequently, the fitted elastic peak was subtracted 
from the experimental data, they were divided by the area of the elastic peak, and the energy 
scale was converted to an energy-loss scale. Finally, the measured spectrum S(T) [in counts 
per channel] was converted to experimental yield y(T) [in reciprocal eV], corresponding to 
Equation ([H]), by division by the channel width AE. Note that due to the dynamical range 
of a typical REELS spectrum measured with good energy resolution, a small misfit of the 
tail of the elastic peak can give erratic excursions in the loss spectra obtained in this way. 
This can be observed in Figure [2] for energy losses below ~5 eV, where a negative excursion 
and a small shoulder right next to it are seen that are due to elimination of the elastic 
peak. The only way to cure this problem is to conduct the experiment with a better energy 
resolution, implying that the primary beam must be monochromatized. 



5. RESULTS. 



The DIIMFP and DSEP retrieved from the spectra displayed in Figure [2] with the second 
order SW-algorithm are presented in Figure [3] as open (DIIMFP) and filled circles (DSEP). 
The solid lines represent the best fit of the data to theory. It is seen that the DIIMFP can 
be perfectly fitted by the employed theory, while the corresponding DSEP agrees reasonably 
with the experimental data, but significant deviations are nonetheless observed. This is 
believed to be attributable to the simplifications concerning the depth dependence of the 
surface excitation process made in the employed theory for surface excitations 18j. The 
Drude-Lorentz parameters giving the best fit between experiment and theory are given 
in Table 1. The binding energy of the most loosely bound core electrons 26| are also 
indicated there and are in good agreement with the energies of the ionization edges observed 
in Figure [31 

A comparison of the real and imaginary part of the dielectric function derived from the 
REELS measurements with the data given in Palik's book [l, 2] is presented in Figure H] and 
[5j Reasonable agreement between these two data sets is observed for all cases for energies 
>5 eV. The error in the present data can become excessively large below 5 eV due to 
problems with the elimination of the elastic peak from the spectra. The error bars in these 
graphs are obtained [8| by assuming that the retrieved DIIMFP predominantly determines 
E2 implying that the uncertainty in £2 is of the order of the uncertainty in the retrieved 
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DIIMFP. The rules of error propagation are then used to estimate the uncertainty in S\. 
As can be seen the uncertainty in e\ is rather large below 20 eV. This is a fundamental 
characteristic of the derivation of optical data from absorption measurements, that mainly 
sample e-i- Within the estimated uncertainty, the two data sets agree satisfactorily, both for 
the real as well as the imaginary part of e. 

The surface (upper panels) and bulk (lower panels) loss functions of Palik's data and the 
present results are compared in Figure [61 For Fe and Ni, the data for the bulk loss function 
used by Tanuma, Powell and Penn 24J is also shown for comparison. The general trend 
observed in these results is that for energies above ~15 eV, the present data show more 
detailed structure in the loss function, while for lower energies, the earlier data seem to be 
more realistic. This is again attributable to the limited energy resolution used in the present 
study and the resulting problem with the elimination of the elastic peak. 

To subject the present optical data to the usual sum-rule checks, the bulk loss function 
was extended above 80 eV by Palik's data. The results for the perfect screening (or "ps") sum 
rule and the Thomas- Reiche-Kuhn (or "f") sum rule are compared with the corresponding 
results based on Palik's and Tanuma's loss functions in Table 2. The ps-sum rule seems to 
be the most important for the present study, since it emphasizes low energies, while the main 
contribution to the f-sum rule comes from the core electrons which are not fully included 
in the present measurements. Except for Ni, the ps-sum rule check for the REELS data is 
in better agreement with the expected value of 1.000 than for the two other sets of optical 
data. The f-sum rule shows deviations of the expected value of atomic electrons of the same 
order of magnitude for all data sets. 

Figure [7] shows the IMFP for the studied materials over the energy range between 50 and 
5000 eV as thick solid line. For comparison, results using the other two sets of optical data 
are also shown as dashed (Palik) and chain dashed (Tanuma) curves. For the calculation of 
the IMFP of all three sets of optical data the software employed by the authors of Ref . 24 ] 
was used [271]. The results predicted by the TPP-2M- formula are indicated by the dotted 
line. The open circles represent the results of elastic peak electron spectroscopy (EPES) 



measurements reported earlier 



25|, 



281 ] . The deviations of the IMFP values based on Palik's 



data set and the present ones are most prominent for Co, which is caused by the lower value 
of Palik's data for the loss function in the region between 20 and 50 eV. For Fe and Ni, 
the mutual agreement between the different IMFP values is satisfactory, except for energies 
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below 100 eV, where the semiempirical TPP-2M formula predicts values for Co and Ni that 
are slightly lower than the other results. The EPES data differ from the other data sets in 
that the experimental elastic peak intensities are interpreted using only the elastic scattering 



cross section as input to the evaluation procedure [25|, [29|], while the other calculations all are 
based on optical data and dielectric response theory. The two approaches are thus in a way 
fundamentally different. Nonetheless, the agreement between the EPES data and the IMFP 
values derived from the dielectric function is reasonable, at least within the experimental 
scatter of the EPES data. 

As a final result, Figure M shows the surface excitation probability extracted from the 
REELS spectra as open circles. This quantity was determined by using the present optical 
data to calculate the normalized DIIMFP and DSEP and by fitting the experimental spec- 
trum to theory via Eqn. ([6]) using (n s ) and the bulk partial intensities as fit parameters. 
Examples of such fits are shown in Figure [2] as thick solid lines. These fits are gene rally 
better than in previous work, where Palik's optica data were used for the same purpose [33] . 
The contribution of electrons that have experienced one bulk, one surface and two surface 
excitations is also indicated in these figures. The solid line in Figure M is a fit of the data 
for the surface excitation probability to Eqn. (j3J), the dotted line represents a fit to another 
functional from of the surface excitation probability that is commonly used 

(n.) = = + = , (18) 

a*vE/ij + l fl*V% + l 



the dashed and chain-dashed line are results by Chen 



311 ] and Kwei et al 32] respectively. 



The quality of these fits is also improved compared to earlier work, in that the va 



the surface excitation probability (n s ) exhibit significantly less scatter than earlier 33] . The 



ues of 



resulting values of the surface excitation parameters a s and a* are given in Table 3. The 
present results for Fe are in excellent agreement with the value of a s = 2.51, reported by 
Chen (dashed line in Figure [8]), while they are also in close agreement with the results quoted 
by Kwei and coworkers [32|. The quality of the fit of the data to Eqn. (jSJ) is generally slightly 
better than the fit to Eqn. f[T8"j) . at least for Co and Ni, although it is still difficult on the 
basis of the present data to decide between the two functional forms for the SEP, Eqn. (J5J) 
and Eqn. (fl8"|) . For this purpose, analysis of REELS experiments at higher energies would 
be required. 
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6. SUMMARY AND CONCLUSIONS. 

An earlier proposed procedure 0, Q] for the simultaneous deconvolution of two REELS 
spectra to provide the energy loss distribution in a single surface and volume electronic 
excitation is simplified using a Pade approximation and applied to experimental data of 
polycrystalline Fe, Co and Ni samples. It is shown that a second order rational fraction 
approximation converges better over an energy loss range of about 150 eV than an 8-th 
order polynomial expansion approximation, allowing one to conclude that the second order 
SW-algorithm can be safely used for practical purposes. Analysis of the expansion coeffi- 
cients provide guidelines on the choice of the optimal experimental parameters to derive the 
DIIMFP and DSEP from REELS spectra and show, moreover, that no prior information 
on surface excitations is needed to perform the deconvolution. The retrieved DIIMFP and 
DSEP were fitted to the corresponding theoretical expressions giving the optical data of 
the studied solids in terms of a set of Drude-Lindhard parameters. Agreement between the 
resulting optical data as well as the IMFP derived from them with values based on optical 
data reported earlier is quite good, supporting the validity of the procedure. The values 
for the surface excitation probability retrieved from the data using the optical constants 
derived in this work are believed to be more reliable than the values reported earlier 33] 
and indicate that the energy and angular dependence of the surface excitation probability 
are described by Eqn. (J5J) rather than by Eqn. (fl8|) . The values for the surface excitation 
parameter are in reasonable agreement with theoretical values. 
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8. APPENDIX: FIRST AND SECOND ORDER PADE RATIONAL FRACTION 
EXPANSION COEFFICIENTS. 



The explicit expressions for the first order bulk expansion coefficients in Eqn. in terms 
of the partial intensities a nbna and fi nb n a of two experimental REELS spectra are given by: 



a b w = An/ A 
a b 01 = -aoi/A 



Ko = 1 

& io = -T^l^o/ffo! - otnPoiPio + «02^i + OioiPwPu - a i/?oi/^2o - am ^ 02 ^ 10 I 

A L Poi > 

h m = -rt\ a loPo2 - aioxonoPu + a^fao + a 10 a n /3 i - a 01 a 2 of3 i - a ° 2ai0 ^ 01 X (19) 
where the determinant A is given by 

A = (aiofti - «oi/5io) (20) 
The surface expansion coefficients read: 

«oo = 



a S io = -Ao/A 
a s oi = Qtio/A 



% = 1 



ttl0/3m/520 



0j> 

2 



&10 — ^|tt20/^01 — a ll/^0l/5l0 — Ctlo/5o2/5lO + «02/5io + «io/^oi/5ii — 

& oi = X2 | a ioA)2 - aoi«io/?ii - «o2«io/?io + «oi/^20 + «oianAo 01 20A0 | 

The expressions for the second order expansion coefficients are somewhat too lengthy to 
reproduce here. They can be easily calculated numerically though, by first establishing the 
polynomial expansion coefficients u h pq and u s pq) as described in Refs. fl2. \\\ up to third order 
(p + q < 3). The coefficients in the denominator of the rational fraction expansion are then 
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found by solution of the homogeneous system of equations: 
= 600 

= u 30 + b w u 20 + b 20 u 10 
= u 2 i + b w u u + b 01 u 20 + b u u w + 6 2 o^io 
= U12 + &10«02 + &oi«n + &11U01 + b 02 ui 
= m 03 + b 0l u 02 + 

= u 22 + b w u 12 + b 01 u 21 + b n u u + 620^02 + b 02 u 20 . (22) 
Finally, the enumerator coefficients are determined by the inhomogeneous system: 

000 = uoo = 

OlO = WlO 

001 = Uoi 

an = Mil + b 10 u 01 + 601^10 

020 = u 20 + 610M10 

002 = u 02 + 601^01- (23) 
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Fe 






Co 




Ni 








E(3p 1/2 ) 






E(3p 1/2 ) 




E(3p 3/2 ) 








52.7 eV 






58.9 eV 




66.2 eV 






A, (eV 2 ) 


7* (eV) 


LOi (eV) 


A, (eV 2 ) 


7, (eV) 


w< (eV) 


A t (eV 2 ) 


7, (eV) 


(eV) 


0.39 


7.54 


0.0 


10.89 


8.16 


0.0 


79.78 


8.16 


0.0 


72.41 


0.50 


1.8 


16.52 


0.50 


2.9 


328.65 


16.19 


1.0 


170.71 


5.95 


4.2 


10.89 


0.50 


4.2 


76.99 


15.64 


2.5 


161.04 


9.53 


10.6 


10.89 


0.50 


5.0 


76.99 


17.07 


3.9 


91.21 


9.12 


19.1 


49.02 


2.51 


6.2 


63.58 


7.04 


14.5 


40.94 


9.74 


27.9 


156.05 


8.15 


9.5 


68.19 


8.39 


23.3 


24.17 


4.78 


34.9 


99.00 


13.94 


15.0 


305.75 


39.04 


41.9 


72.99 


3.70 


55.3 


146.69 


13.57 


19.9 


84.22 


14.48 


49.8 


64.36 


8.06 


60.8 


10.89 


3.10 


25.5 


164.52 


29.48 


64.2 


292.44 


31.28 


74.0 


99.00 


26.18 


35.7 


98.01 


6.92 


68.7 








15.39 


5.62 


38.8 


7.59 


1.35 


76.9 








10.89 


5.74 


50.8 














60.53 


3.92 


60.9 














54.20 


8.11 


68.3 









TABLE I: Drude-Lindhard parameters for Fe, Co and Ni determined in the present work 





Fe 


(Z=26) 


Co 


(Z=27) 


Ni 


(Z=28) 


Ref. 


f-sum 


ps-sum 


f-sum 


ps-sum 


f-sum 


ps-sum 


REELS 


24.0 


1.007 


24.2 


0.925 


29.4 


1.030 


Palik 


24.0 


0.943 


24.2 


0.845 


26.7 


1.010 


Tanuma 


23.5 


1.110 






27.3 


1.050 



TABLE II: Comparison of f- and ps-sum rule checks for different sets of optical data. The REELS 
data above 100 eV were extended to 30 keV using the data set of Palik 
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a s (Eqn. ©) 


a*(Eqn. ®) 




(eV l l 2 ) 


(eF- 1 ^) 


Fe 


2.53 


0.34 


Co 


2.79 


0.30 


Ni 


2.84 


0.30 



TABLE III: Surface excitation parameter a s , according to Eqn. (JSJ) and a*, according to Eqn. (|18j) 
determined from the surface excitation probabilities shown in Figure [8] 
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Energy loss T (eV) 
0.03 l 




Energy loss T (eV) 

FIG. 1: Illustration of the successive stages of the deconvolution of REELS spectra of Fe (measured 
at 1000 and 3400 eV) using an eighth order polynomial approximation (a.) and a second and first 
order Pade rational fraction approximation (Eqn. (|8|)) (b). The open circles in (b) are the result 
of the 8-th order polynomial expansion Eqn. ([7]), the solid curves labelled (1) and (2) are the first 
and second order Pade approximation, respectively. The curves shown represent the DIIMFP, the 
corresponding results for the DSEP are not shown for clarity, but behave in a similar way. 
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(c) 
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FIG. 2: Experimental loss spectra for Fe, Co and Ni used in the present work (noisy curves). The 
thick solid curve is a fit of these data to Eqn. ([6]), using the retrieved optical data to describe the 
DIIMFP and DSEP. The thin solid curves labelled " Is" , " lb" and " 2s" represent the contribution of 
electrons having experienced one surface, one bulk and one and two surface excitations respectively. 
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FIG. 3: Normalized differential surface excitation probabililty (DSEP, w s (T)) and normalized 
differential inverse inelastic mean free path (DIIMFP, Wb(T)) retrieved from the experimental data 
by means of Eqn. Q. The solid line is obtained by simultaneously fitting these data to theory 
for a set of Drude-Lorentz fit-parameters for the dielectric function, as described by Eqn. (I16j) . 
The data for the surface excitation probability were scaled by a factor of 0.25 in order to facilitate 
comparison, (a) Fe; (b) Co; (c) Ni. 
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FIG. 4: Real part of the dielectric function, ei, of Fe, Co and Ni retrieved from the REELS data 
(open circles) compared with Palik's data lj, |2| (solid curves), (a) Fe; (b) Co; (c) Ni. 
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FIG. 5: Imaginary part of the dielectric function, ei , of Fe, Co and Ni retrieved from the REELS 
data (open circles) compared with Palik's data [lj, |2( (solid curves), (a) Fe; (b) Co; (c) Ni. 
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FIG. 6: Volume and surface loss function of Fe, Co and Ni retrieved from REELS data (open 

circles) compared with Palik's optical data (dashed curves). For the volume loss function of Fe 

and Ni, Tanuma's data are also shown as solid curves, (a) Fe; (b) Co; (c) Ni. 
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FIG. 7: Inelastic mean free path of Fe, Co and Ni calculated from the optical data shown in Figure[6] 
(solid curves); Open circles represent experimental results measured with elastic peak electron 



spectroscopy (EPES) 



spect; 
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281 ]; Dashed curves represent calculations based on Pali 



The dotted curves represents the values predicted by the TPP-2M formula 



f's optical data 
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FIG. 8: Surface excitation probability as a function of the energy for (a) Fe; (b) Co; and (c) 
Ni (open circles). The solid and dashed lines are fits of these data to Eqn. ([5] and Eqn. (|18p 



3l| which is 



respectively). The dashed line shown for Fe is the theoretical dependence by Chen 
almost impossible to distinguish from the present data. The chain-dashed curves for Fe and Ni 
represent the results by Kwei et al [321 ] . 
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